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Abstract 

We present exact solutions of the incompressible Navier-Stokes equations in a background linear 
shear flow. The method of construction is based on Kelvin's investigations into linearized distur- 
bances in an unbounded Couette flow. We obtain explicit formulae for all three components of 
a Kelvin mode in terms of elementary functions. We then prove that Kelvin modes with paral- 
lel (though time-dependent) wavevectors can be superposed to construct the most general plane 
transverse shearing wave. An explicit solution is given, with any specified initial orientation, profile 
and polarization structure, with either unbounded or shear-periodic boundary conditions. 
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In 1986 Craik and Criminale [l| presented a class of exact solutions of the Navier-Stokes 
equations which were wavelike disturbances in background shear flows. Since then these 
solutions have proved extremely useful in the study of astrophysical and atmospheric fluid 
dynamics; a very useful collection of exact solutions can be found in [2j. The approach 
taken in |l[ was a generalization of a century-old method invented by Kelvin |3j to study 
linearized perturbations of Couette flows. These shearing wave solutions, also referred to as 
Kelvin modes, have time-dependent wavevectors and amplitudes. Although a single Kelvin 
mode is an exact solution of the full Navier-Stokes (NS) equations, it has been remarked 

n 

PJ that until about 1965 there seems to be no evidence that this was so recognized; in fact, 
the first published mention is as late as 1983 4]. Moreover, an explicit formula has been 

nfi 

published |l|, [3] for only one of the three components of the disturbance. 

In this Letter we present exact solutions for all three components of the velocity field of 
a Kelvin mode, in closed form using only elementary mathematical functions. We identify 
a subset of these modes whose wavevectors — though time-dependent — remain parallel 
to each other for all time. These are used to synthesize the most general plane transverse 
shearing wave, which can have any specified initial orientation, profile and polarization 
structure, with either unbounded or shear-periodic boundary conditions. 

Let (ei, e 2 , £3) be the unit basis vectors of a Cartesian coordinate system in the laboratory 
frame. Using notation x = (xi,x 2 ,x^) for the position vector and t for time, we write the 
total fluid velocity as (Sx\e 2 + v), where S is the rate of shear parameter and v(x, t) is the 
incompressible disturbance (V-v = 0) which obeys the NS equations: 

(d t + Sx x d 2 )v + Svie 2 + (v-V)v = -Vp + vV 2 v , 

V 2 p = -V-[(«-V)«] - 2Sd 2 Vi. (1) 
We seek a solution in the form of a single Kelvin mode 

v k {x } t) = Re{A(fc,t) exp [i k sh (t)-x] } , 

p k (x,t) = Re{ip(k,t) exp [ik sh {t)-x]} , (2) 
where the time-dependent sheared wavevector, k sh (t), has components 
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with k = (fci, fc 2 , £3) being a constant wavevector. Our task now is to determine the ampli- 
tudes A(fe, t). Incompressibility requires that k sh (t)-A(k, t) = 0. Therefore, when eqns. (J2]) 
and (JSJ) are substituted in eqns. (CQ), the nonlinear term, (v«V) t> vanishes because 



(A-V) exp [ifc sh (t)-;c] = (ifc sh (t)-A) exp [ifc sh (t)-a;] = 0. 

The pressure can be eliminated by using the second of eqns. 
where |fc sh (t)| 2 = [(h - Stk 2 f + k% + kf\ . Then A satisfies 
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k 2 k sh (t) 



d t A + SA& = 2S 



, . A x - v\k sh (t)\ 2 A. 
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We now obtain explicit solutions for A. To do this, define a new amplitude variable, a(fc, t), 
by 



A(fc,t) = G v (k,t) a{k,t) . 
where G u (k,t) is a Fourier-space viscous Green's function, 
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G u (k,t) = exp 



-v I ds \k sh (s)\ 2 
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When eqns. and ([7]) substituted in eqn. ([5]), we obtain the following equations for the 
three components of a(k,t) : 



d t a x - 2S 
d t a 2 - 2S 



- Stk 2 ) k 2 
(k ± - Stk 2 ) 2 + k 2 + k 2 
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d t a 3 - 2S 



(h - Stk 2 y + k 2 + k 2 
k 2 h 



(ki - Stk 2 y + k 2 + k 2 



ai = . 



ai = . 



ai = 



(8) 
(9) 
(10) 



Eqn. OH]) can be solved to get an explicit expression for ax(k,t): 



3 



at(k,t) 



k 2 

'/.-,)'- + k\ 



k\ 



ai(fe,0) 



(11) 



{k\ — OLK2) 

which expression is given in |3j. When this is substituted in eqns. OH]) and fllOl) . the latter 
can be integrated to obtain expressions for a 2 (k,t) and a 3 (k,t). However, neither Kelvin 
nor anyone else, to the best of our knowledge, have published explicit formulae for these two 
components. Thus we were pleasantly surprised to find that a 2 (k,t) and as(k,t) could be 
expressed entirely in terms of elementary functions: 



a 2 (k,t) = a 2 (k,0) 
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a 3 (k,t) = a 3 (k,0) - 
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Incompressibility requires that k sh (t)-a(k,t) = 0, which is guaranteed if the initial condi- 
tions are chosen such that k-a{k, 0) = . When eqns. ©, (J7J), (TTTj) and ([121) are substituted 
in eqn. (T5]), we obtain the full velocity field of a single Kelvin mode. Until now we have 
considered an unbounded flow. However, in numerical simulations of the local dynamics 

n n 

of differentially rotating discs in astrophysical systems [5|, |6[, it is customary to employ 
"shear-periodic" boundary conditions. Let us define sheared coordinates by 



x 



sli 



Xi , 



x" 2 = x 2 — Stx\ , xf 1 = £3 . (13) 

These may be thought of as the Lagrangian coordinates of fluid elements that are carried 
along by the background shear flow. A function is said to be shear-periodic when it is a 
periodic function of (xf 1 , x| h , xf 1 ) with periodicities (L 1; L 2 , L 3 ), respectively. The phase of 
the function can be written as k sh (t)-x = k-x sh . Therefore, a shear-periodic Kelvin 
mode has wavevectors k £ {2-nmxj L\ } 27nn 2 /L 2 , 27rm3/L 3 ), where the m 8 take any integer 
values. 
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We now use the explicit expressions obtained for the Kelvin modes to construct the most 
general plane transverse shearing wave. Let us consider two Kelvin modes, Vu(x,t) and 
Vu'(x,t) corresponding to wavevectors k and k', which are parallel to each other but could 
differ in magnitudes. Using eqns. Q, we see that the corresponding sheared wavevectors, 
k sh (t) and k sh (t), are also parallel to each other for all time. Incompressibility implies 
that Vfc(x,t) and v^i{x,t) are perpendicular to k sh (t) and k sh (t) for all time. So, if we 
superpose vj^(x,t) and Vu'(x,t), the nonlinear term in the NS equations vanishes, because 
the superposed velocity field remains parallel to the wavefronts. Thus the superposition of 
an arbitrary number of Kelvin modes, all with wavevectors parallel to each other, is an exact 
solution of the NS equations. 

Let us choose a unit vector h = (ni, n 2 , n s ), and define the sheared (non-unit) vector 
n sh (t) by 

< = (m - Stn 2 ) , nf = n 2 , nf = n 3 . (14) 

Superposing all Kelvin modes with wavevectors q = qn, where — oo < q < oo, we obtain 
an exact plane-wave solution of the NS equations with wavefronts perpendicular to n sh (t): 



Vi(x,t) 



dq 
2^ 



G v (qh,t)Wi(q) exp [iqn sh (t)-x] + 
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Fi (n sh (t)) - Fiih) 



p- G u (qh,t) WM exp [iqn s \t)-x] , (15) 



where the dimensionless and scale-invariant functions, Fi(Q), are defined by 
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For shear-periodic boundary conditions, the integral over q in eqn. (Tl5|) should be replaced 
by an appropriate sum. The W(q) are Fourier-space initial conditions corresponding to the 
a(k, 0) of eqns. ( flTl) and (fl2l) . and must satisfy the incompressibility condition, n-W(q) = 
. They are determined by the initial profile and polarization stucture of the plane wave. 
At t = 0, the wavefronts are perpendicular to n, so we write v(x,0) = W(n-x), where 
h-W = 0. Note that the only constraint on the initial condition, W, is that it is a 
vector field that is perpendicular everywhere to the unit vector n; otherwise it is a quite 



arbitrary function of its one argument. Thus, no restriction need be placed on the initial 
profile and polarization structure of the initial conditions. Given W(y), we can determine 
W(q) = dy W(y) exp [— iqy] , and use this in eqn. (JT5|) to calculate v(x, t). 

Eqn. (|T5|) is a complete solution for a general plane shearing wave, expressed in terms of 
a Fourier integral. However, it is physically more transparent to rewrite the right side in 
terms of real-space quantities. To do this, we must introduce the real-space viscous Green's 
function, whose natural definition is with respect to the sheared coordinates [TJ: 

/c\k ~ 
— G u (k,t)exp[ik.x sh }. (17) 

The properties of this function are discussed in S], where it is shown that it takes the 
form of a sheared heat kernel, which is an anisotropic Gaussian function of x sh with time- 
dependent coefficients; all the principal axes increase without bound and rotate against the 
direction of the background shear. Noting that n sh (t)-x = n-x sh , we can write the general 
form of the plane shearing wave as 



Vi (x,t) = I d a £G v (£,t)W i (n-[x ai (t)-€]) + 

J d^G v (tt)W 1 (n.[x sh (t)-$}) (18) 

As an illustrative example let us consider the following initial condition, corresponding to 
a polarized wavepacket with wavevector pointing along the x 2 -a.xis: n = e 2 , Wi(x 2 ) = 
W exp [-xl/2a 2 ] smkx 2 , W 2 = , W 3 (x 2 ) = hW exp [-xl/2cr 2 ] cos kx 2 , where -1 < h < 
1 . The wavepacket is linearly polarized when h = 0, and right/left circularly polarized when 
h = ±1; other values of h correspond to different degrees of elliptical polarizations. At a 
later time, the wavevector has components nf 1 = —St , n| h = 1 , nf 1 = . Since both Wi and 
G u (£, t) are Gaussian functions, the integrals in eqn. (1181) can be performed analytically and 
v(x,t) evaluated explicitly. We present the results graphically in Fig.([T]) for two cases, one 
linearly polarized and the other right circularly polarized. As the wavepackets are sheared, 
they undergo transient amplification due to the combined action of shear and viscosity, and 
at late times suffer viscous damping. 
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In conclusion, we have constructed exact solutions of the Navier-Stokes equation with a 
background linear shear flow. All three components of the velocity field of the Kelvin modes 
are given in closed form using only elementary mathematical functions. It is demonstrated 
that, when Kelvin modes with parallel wavevectors are superposed, they remain exact solu- 
tions. We give in explicit form the most general plane transverse shearing waves, with any 
specified initial orientation, profile and polarization structure, with either unbounded or 
shear-periodic boundary conditions. These represent the local structure of any disturbance 
in general shear flows, and can therefore be expected to find many applications in the theory 
and simulations of astrophysical and atmospheric flows. 
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FIG. 1. (color online) Evolution of plane transverse shearing wavepackets. The polarization struc- 
ture of the velocity field is indicated on some sections of the plane wavefronts. The parameters 
values used are S = — 1, v = 1, Wo = 1, a = 10 and k = 1. (a) and (b) Linearly polarized (h = 0) 
wavepackets at times t = and t = 1. (c) and (d) Right circularly polarized (h = 1) at times t = 
and i = 1. 
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